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We present a study of mechanical, electronic and magnetic properties of two dimensional (2D), 
monolayer of silicon carbide (SiC) in honeycomb structure and its quasi ID armchair nanoribbons 
using first-principles plane wave method. In order to reveal dimensionality effects, a brief study of 
3D bulk and ID atomic chain of SiC are also included. Calculated bond-lengths, cohesive energies, 
charge transfers and band gaps display a clear dimensionality effect. The stability analysis based on 
the calculation of phonon frequencies indicates that 2D SiC monolayer is stable in planar geometry. 
We found that 2D SiC monolayer in honeycomb structure and its bare and hydrogen passivated 
nanoribbons are ionic, non-magnetic, wide band gap semiconductors. The band gap is further 
increased upon self-energy corrections. The mechanical properties are investigated using the strain 
energy calculations. The effect of various vacancy defects, adatoms and substitutional impurities on 
electronic and magnetic properties in 2D SiC monolayer and in its armchair nanoribbons are also 
investigated. Some of these vacancy defects and impurities, which are found to influence physical 
properties and attain magnetic moments, can be used to functionalize SiC honeycomb structures. 

PACS numbers: 73.22.-f, 75.75.+a, 63.22.-m 



I. INTRODUCTION 

Owing to its exceptional thermal and physical 
properties,^ silicon carbide (SiC) is a material, which is 
convenient for high temperature and high power device 
applications. Because of its wide band gap, SiC bulk 
structure has been a subject of active study in optical 
and optoelectronic research. Unlike the polymorphs of 
carbon, SiC is a polar material. In spite of the fact that 
both constituents of SiC are Group IV elements, charge 
is transferred from Si to C due to higher electronegativity 
of C relative to Si atom. 

Bulk SiC has six commonly used stacking configura- 
tions denoted as 3C (zincblende), 2H (wurtzite), 4H, 6H, 
15R and 21R. Lubinsky et reported optical data re- 
lated with indirect transitions, dielectric function and re- 
flectivity of 3C SiC using first-principles Hartree-Fock- 
Slater method. A more comprehensive study^ using OL- 
CAO (orthogonalized linear combination of atomic or- 
bitals) method comprises the calculations of lattice con- 
stants, electronic band structure and optical properties 
of all six stacking configuration of SiC. 

As for SiC in lower dimensionality, Si02 coated SiC 
nanowires^ were synthesized and showed favorable pho- 
tocatalytic behavior. A theoretical work on hydrogen 
passivated SiC nanowires 5 provided the energy bands 
both using local density approximation within Density 
Functional Theory (DFT) and sp 3 s* LCAO tight bind- 
ing (TB) methods. SiC-ZnS core-shell structures were 
also fabricated. e Zincblende SiC nanoparticles were syn- 
thesized by carbothermal reduction methodP' Band gap 
of zincblende nanoparticles were estimated to be around 
3 eV from photoluminesence measurements. With a sim- 
ilar carbothermal method, microribbons^ with widths in 
the range of 500 nm - 5 /im and thickness of 50-500 nm 



were synthesized. 

SiC is freq uently used as a substrate to grow other 
materials jr ew layers of graphene was also grown on 
SiCPSiC clusters (Si„C„, n=l-10) were investigated^ 
using DFT. With the aim of developing a material 
for future nanoelectronic applications, binding energy, 
HOMO-LUMO gap, Mulliken charge, vibrational spec- 
trum and ionization potential of Si„C n clusters are re- 
vealed. 

Earlier, planar honeycomb structure of graphite was 
exfoliated and its physical properties were analyzed.^KS 
While graphene is a strictly planar crystal, the planar 
honeycomb structure of Si is unstable, but it is stabilized 
through puckeringfED Since the honeycomb structure is 
common to both C and Si, one expects that stable 2D 
SiC in honeycomb structure can be synthesized. 

In this paper, a comprehensive analysis of the atomic, 
electronic and magnetic properties of 2D monolayer of 
SiC honeycomb structure and its bare and hydrogen pas- 
sivated armchair nanoribbons (A-SiCNR) are carried out 
using first-principles calculations. In spite of the fact 
that 2D SiC monolayer is not synthesized yet, this study 
demonstrates its stability based on reliable theoretical 
methods. Furthermore, various mechanical, electronic 
and magnetic properties are revealed. We started with 
the discussion of 3D zincblende and wurtzite crystals, 
as well as SiC atomic chain as an ultimate ID system; 
we presented an analysis of optimized atomic structures 
with corresponding phonon dispersion curves and elec- 
tronic energy band structures and effective charges. Then 
we provided an extensive analysis of 2D and quasi- ID 
(nanoribbon) SiC in terms of the optimized atomic struc- 
tures and their stability, electronic and magnetic struc- 
tures. We revealed elastic constants, such as in-plane 
stiffness and Poisson's ratio. Having obtained the results 



for ID, 2D and 3D structures, we presented a compre- 
hensive discussion of dimensionality effects. Then we in- 
vestigated the effect of vacancy defects (such as Si and C 
vacancy, Si+C vacancy and C-Si antisite defect) on the 
electronic and magnetic properties of single layer SiC and 
its armchair nanoribbons. Furthermore, we showed that 
SiC can be functionalizcd through adsorption of a foreign 
atom to the surface of 2D SiC or through substitution of 
either C or Si with a foreign atom. It is revealed that 2D 
SiC and its ribbons provide unusual physical properties, 
which are absent in 3D SiC crystals. For example, while 
various allotropic forms of SiC including its honeycomb 
structures are normally non-magnetic semiconductors, a 
Si vacancy gives rise to spin polarization. Significant vari- 
ation of the band gap of narrow A-SiCNR's with their 
widths may be crucial in designing optoelectronic nan- 
ode vices. 



II. MODEL AND METHODOLOGY 

We have performed first-principles plane wave calcula- 
tions within DFT using PAW potentials!-^ The exchange 
correlation potential has been approximated by Gener- 
alized Gradient Approximation (GGA) using 
functional both for spin-polarized and spin-unpolarized 
cases. For the sake of comparison, the calculations are 
also carried out using different potentials and exchange- 
correlation approximations. All structures have been 
treated within the supercell geometry using the periodic 
boundary conditions. A plane- wave basis set with kinetic 
energy cutoff of 500 eV has been used. A vacuum spacing 
of 12 A hinders the interactions between SiC monolay- 
ers in adjacent supercells. In the self-consistent structure 
optimizations, the Brillouin zone (BZ) is sampled by, re- 
spectively (5x5x5), (llxllxl) and (llxlxl) special k- 
points for 3D bulk, 2D honeycomb and ID nanoribbons 
of SiC. Further relaxation is made with (llxllxll), 
(31x31x1) and (25x1x1) special k-points in order to 
find the final structure. All atomic positions and lattice 
constants are optimized by using the conjugate gradi- 
ent method, where the total energy and atomic forces 
are minimized. The convergence for energy is chosen as 
10~ 5 eV between two steps and the maximum Hellmann- 
Feynman forces acting on each atom is less than 0.04 
eV/A upon ionic relaxation. The pseudopotentials cor- 
responding to 4 valence electrons of Si (Si :3s 2 3p 2 ) and 
C (C:2s 2 2p 2 ) are used. Numerical plane wave calcu- 
lations are performed by using 

VASPPMI p ar t f the 
calculations have also been repeated by using SIESTA 21 
software. The cohesive energy of any SiC structure is 
found as E c = E T [SiC] - Et [Si] - E T [C] in terms of 
the optimized total energy of any SiC structure, and the 
spin-polarized total energies of free Si and C atoms, all 
calculated in the same supercell using the same parame- 
ters. Phonon calculations were carried out using PHON 
program 2 ^ implementing force constant method. GWq 
calculations^ are again handled by VASP. 
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Figure 1: (Color online) Atomic structures, electronic energy 
bands and dispersion of phonon modes of linear and wide 
angle zigzag SiC atomic chains. Ec and Eg are cohesive and 
band gap energies, respectively. Si and C atoms are shown by 
blue/large and brown/small balls, respectively. Zero of the 
energy is set to the Fermi energy, Ep. 



III. ID ATOMIC CHAIN AND 3D BULK 
CRYSTAL OF SIC 

In this section, we present a brief discussion of ID SiC 
atomic chain and 3D bulk crystal based on our struc- 
ture optimized total energy and phonon spectrum cal- 
culations. Studies on SiC bulk lattice and atomic chains 
already exist in the literature! 2 * 3 * 24 ! However, our purpose 
is to carry out calculations with same parameters as used 
in 2D monolayer SiC honeycomb structure and provide 
a consistent comparison of dimensionality effects. 



A. ID SiC Chains 

Earlier, the first theoretical study of atomic chains of 
Group IV and III-V binary compounds were reported 
by Senger et al¥^ They examined SiC atomic chain as 
a function of lattice parameter and found that the wide 
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Table I: Si-C bond length, d; lattice constant, a; bond angle, 
a; charge transfer from Si to C, 8q; band gap, Eg; and cohe- 
sive energy, Ec values for two different types of SiC chains. 



Type 


d (A) 


a (A) 


a (degree) 8q (e) E G (eV) E c (eV) 


Linear 


1.649 


3.298 


180 2.28 2.00 8.962 


Wide Zigzag 


1.673 


3.268 


155.2 2.15 1.18 8.968 



zigzag atomic chain of SiC with bond angle of ^147° 
is energetically more favorable than the linear and nar- 
row angle zigzag chains. Present calculations find that 
the atomic chains of SiC are non-magnetic. Calculated 
structural parameters, cohesive energies, band gap and 
phonon modes of linear and zigzag atomic chains, which 
are relevant for the present study are given in Table [I] 
The charge transfer from Si to C is calculated to be <5q= 
2.28 electrons using the Bader analysis! 2 ^* Phonon modes 
calculated with force constant method have imaginary 



frequencies. In Fig. [T] two acoustical and two optical 
branches of linear SiC chain have imaginary frequencies. 
Also wide angle zigzag SiC chain has one optical and one 
acoustical branches with imaginary frequencies. These 
results indicate that free standing SiC chains are not sta- 
ble. We note that carbon and BN ato mic chain s are found 
to be stable and have linear structure! 24 * 26 * 27 ! Stability of 
linear chain structure is assured by 7r-bonding between 
adjacent atoms. 



B. 3D SiC Crystals 

Our work on bulk SiC includes wurtzite (wz) and 
zincblende (zb) structures. Atoms in wz- and zb-SiC 
are four fold coordinated through tetrahedrally directed 
suborbitals. Calculated structural parameters, cohesive 
energies, energy band structures and phonon modes are 
given in Fig. [2] Zincblende SiC structure in symmetry 
has cubic lattice constants, ai = a-i = = 3.096 A. Si-C 
bond distance d is 1.896 A. Each Si (C) is connected to its 
four nearest neighbor C (Si); four Si-C bonds are equal. 
Charge transfer from Si to C is 5q= 2.59 electrons calcu- 
lated via Bader analysis.^ While the GGA band gap is 
1.41 eV, it increases to 2.40 eV after GWq corrections. As 
for wz-SiC crystal, the hexagonal lattice constants of the 
optimized structure in equilibrium are a\ = a 2 = 3.091 
A, c/a — 1.642. The small deviation of c/a from the ideal 
value of 1.633 imposes a slight anisotropy on the lengths 
of tetrahedrally directed Si-C bonds. While the length of 
three short bonds is 1.893 A, the fourth bond is slightly 
longer and has the length of 1.907 A. Charge transfer 
from Si to C is Sq= 2.63 electrons. The GGA band gap 



is 2.36 eV, but it increases to 3.32 eV after GWq correc- 
tion. The calculated structural parameters and energy 
band gaps are in reasonable agreement with the earlier 
calculations and experimental measurementsP^H I n par- 
ticular, the band gap values of 2.39 and 3.33 eV for zb 
and wz SiC, respectively are in excellent agreement with 
the present GWq corrected values. 

The frequencies of phonon modes and their disper- 
sions are calculated for zb and wz crystals by direct (or 
force constant) method.^ At long wavelengths near the 
r point, the electric field due to dipoles is critical for 
phonon modes. This effect lifts the degeneracy between 
longitudinal and transverse optical modes. However the 
splitting (known as LO-TO splitting) cannot be observed 
with the direct method used in the present study. There- 
fore, in Fig. [2] the highest and second highest optical 
branches become degenerate at V point. Present result s 
are in agreement with earlier phonon calculations ! 29 * 3 ^ 

IV. 2D SIC HONEYCOMB STRUCTURE 



Two-dimensional monolayer of SiC hexagonal structure 
of SiC is optimized using periodically repeating supercell 
having 12 A spacing between SiC planes. The minimum 
of total energy occurred when Si and C atoms are placed 
in the same plane forming a honeycomb structure. The 
magnitude of the Bravais lattice vectors of the hexago- 



nal lattice is found to be a\ = a 2 = 3.094 A (see Fig. I), 
and the Si-C bond length to be d =1.786 A. The pla- 
nar structure of 2D SiC is tested by displacing Si and C 
atoms arbitrarily from their equilibrium positions by 0.5 
A and then reoptimizing the structure. Upon optimiza- 
tion, the displaced atoms returned to their original posi- 
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Figure 2: (Color online) Optimized atomic structure with 
relevant structural parameters, corresponding energy band 
structure and frequencies of phonon modes of 3D bulk SiC in 
zincblende and wurtzite structures. Zero of energy of the band 
structure is set at the Fermi level, and band gap is shaded. 



Table II: Si-C bond length, d; lattice constant, a; band gap, 
Eg; band gap corrected by GWo, Egw ; cohesive energy, 
Ec values for 2D monolayer of SiC in honeycomb structure 
calculated with different potentials 



Potential 


d (A) 


a (A) 


E G 


E G (+GW ) 


E c (eV) 


PAW+GGA 


1.786 


3.094 


2.530 


3.90 


11.944 


PAW+LDA 


1.770 


3.070 


2.510 




13.542 


US+GGA 


1.776 


3.079 


2.542 




11.973 


US+LDA 


1.759 


3.048 


2.532 




13.471 



tions in the same plane implying the stability of planar 
structure. Further details on the stability of this struc- 
ture will be provided with phonon calculations at the end 
of this section. Two dimensional monolayer SiC is found 
to be a semiconductor with a band gap of 2.53 eV within 
GGA. Furthermore, in Table [TTJ we present lattice con- 
stant, bond length, cohesive energy, energy gap values 
of 2D SiC monolayer calculated with different potentials. 
Since DFT usually underestimates the band gap of semi- 
conductors, we also corrected the GGA band gap using 
GWo correction and found it to be 3.90 eV. The charge 



transfer from Si to C in 2D SiC is calculated to be <5q= 
2.53 electrons. The Si-C bond length and the band gap 
values in the first row obtained by using GGA+PAW are 
in agreement with earlier DFT calculations.^ 

In addition to structural and electronic properties, 
we also investigated the mechanical properties of 2D 
SiCJSmHU One can specify the mechanical properties of 
SiC honeycomb sheet in terms of uniaxial strain, e = 
Ac/c, c being the lattice constant; Poisson's ratio v — 

etrans/eaxial and in-plane stiffness, C= ^(%JH- Here 
Aq is the equilibrium area of the system and E s is the 
strain energy calculated by subtracting the total energy 
of the strained system from the equilibrium total energy. 
To calculate the elastic constants of monolayer SiC hon- 
eycomb structure we switch to rectangular unit cell and 
consider a large (8x4) supercell comprising 32 primitive 
unit cells. In the harmonic region E s (e) is first calcu- 
lated on a 2D grid. The numerical data is then fitted 
to the expression, E s = Aie 2 + ^be 2 , + A^€ x € y \ where e x 
and e y are the small strains along x- and y-directions. In 
the harmonic region, as a result of the isotropy A\ = A?,. 
More details concerning the calculation of C and v can be 
found in references ! 33 * 34 ! The calculated in-plane stiffness 
of SiC honeycomb structure is found to be 166 (J/m 2 ). 
This is almost half of the in-plane stiffness of graphene 
(namely 335 J/m 2 ), but more than twice the in-plane 
stiffness of silicene (62 J/m 2 ). Also the Poisson's ratio 
of SiC is calculated to be 0.29. Consequently, 2D SiC 
monolayer is a stiff material, but less stiff than graphene 
and BN having similar honeycomb structures. 

Similar to ID and 3D SiC, frequencies of phonon modes 
of and the dispersions fi(k) of 2D monolayer of SiC in 
planar geometry are calculated using direct method.^ 
Forces were found by displacing a single atom in a 7x 7x 1 
supercell. We use a small displacement in order to stay 
in the harmonic region. We increased default grids used 
by VASP until calculations converge. The lowest acous- 
tical mode, which is called as the out of plane, ZA mode 
is vulnerable to instability. In rough meshes this mode 
gets imaginary frequencies near the T-point, but it can 
be overcome by refining the mesh along the z-axis (per- 
pendicular to plane) as much as possible. This way, force 
in that direction is calculated more rigorously. Since the 
frequencies of all modes are positive in BZ, it is concluded 
that planar, 2D SiC monolayer in honeycomb structure 
is stable. 

A. Dimensionality effects 



In Table |III[ we compare the variation of the effective 
charge on Si and C atoms, namely Z* Si and Z* c respec- 
tively; charge transfer from Si to C, Sq = 4 — Z* Si \ Si-C 
bond length d; lattice constant; energy band gap, GWo 
corrected band gap and cohesive energy per Si-C, cal- 
culated for SiC in different dimensionalities. It should 
be noted that the length of Si-C bonds of 2D SiC hon- 
eycomb structure is smaller than that in the 3D bulk 
(wz, zb) crystals, but larger than that in zigzag atomic 
chains. Here we see that the dimensionality effect is re- 
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Table III: Bonding types; Si-C distances, d; lattice constants, a; charge transfers, Sq; band gaps, Eg with GWo corrections; 
and cohesive energies, E c for comparison of SiC polymorphs 



structure 


Bonding 


d (A) 


a (A) 


Sq (e) 


Zh 


Z*c 


E G (eV) 


E G (+GWo) 


E c (eV) 


Linear Chain 


sp 


1.649 


3.298 


2.28 


1.72 


6.28 


2.00 




8.923 


Wide Zigzag Chain 


sp 


1.673 


3.268 


2.15 


1.85 


6.15 


1.18 




8.963 


2D Honeycomb 


sp 2 


1.786 


3.094 


2.53 


1.47 


6.53 


2.53 


3.90 


11.940 


Zincblende 


sp 3 


1.896 


3.096 


2.59 


1.41 


6.59 


1.41 


2.40 


12.939 


Wurtzite 


sp 3 


1.893 (3), 1.907 (1) 


3.091 


2.63 


1.37 


6.63 


2.36 


3.32 


12.933 



(a) Structure (b) Energy Bands (c) Phonon Dispersions 




Figure 3: (Color online) (a) atomic structure, (b) energy band structure and (c) phonon modes of 2D SiC in honeycomb 
structure. Large/blue and small/brown balls indicate Si and C atoms respectively. The primitive unit cell is delineated. The 
zero of energy in the band structure is set to the Fermi level. 



fleeted to the strength of the bonding through sp n hy- 
bridization, where n coincides with the dimensionality. 
While sp 2 hybrid orbitals of 2D planar honeycomb struc- 
ture form stronger bonds than tetrahedrally coordinated 
sp 3 orbitals of 3D bulk, they are relatively weaker than 
sp hybrid orbitals of ID chain. Accordingly, d is shortest 
in ID chain, and longest in 3D zb structures, and is in- 
termediate in 2D monolayer. The cohesive energy Eq in- 



V. BARE AND HYDROGEN PASSIVATED SIC 
NANORIBBONS 



In this section, we consider bare and hydrogen pas- 
sivated armchair SiC nanoribbons. These nanoribbons 
are specified according to their widths specified in terms 
of N number of Si-C basis in their unit cells. Hence, 
A-SiCNR(TV) indicates armchair SiC nanoribbons hav- 
ing TV Si-C pairs in their unit cell. We have analyzed 
A-SiCNR (both bare and H-passivated) from iV=5 to 
21. A-SiCNR's with odd numbered N have reflection 
symmetry with respect to their axis. Bare armchair SiC 
nanoribbons are ferromagnetic in ideal honeycomb form. 
However, upon structural relaxation, reconstruction oc- 
curs at the edges resulting in a considerable gain of energy 
and the structure becomes non-magnetic. 



creases with dimensionality, since the number of nearest 
neighbors increases. Effective charge or charge transfer 
between cation and anion also varies with dimensionality. 
For example, the charge transfer calculated with Bader 
analysis increases with increasing dimensionality. While 
the energy band gap Eq does not show a regular trend 
with dimensionality, the band gap of 2D monolayer is 
wider than those of 3D crystals. 

I 

Band gaps of all A-SiCNR(N) increase by ~ 0.7 eV 
upon H-saturation of dangling bonds. The Si-H and C-H 
bonds formed after hydrogenation have lengths of 1.49 A 
and 1.09 A, respectively. The effective charges on Si and 
C edge atoms changes after H-passivation. While Si-C 
bonds at the edges are shorter than the Si-C bonds in 
2D SiC monolayer by 0.09 A, upon H-saturation these 
bonds get slightly longer, but are still shorter than the 
regular Si-C bonds by 0.05 A. Hence, reconstruction of 
atomic structure at the edge exist in each case, but is 
more pronounced in the bare nanoribbons. 

Here we consider A-SiCNR(9) as a prototype and ex- 
amine its band structure. The bare A-SiCNR(9) is an 
indirect band gap semiconductor. Two bands at the con- 
duction band edge are composed of edge states, which are 
split due to edge-edge interactions. These bands are re- 
moved upon H-saturation of dangling bonds of atoms at 
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Figure 4: (Color online) Energy band structure of the bare 
(a) and hydrogen saturated (b) armchair SiC nanoribbons, 
A-SiCNR(iV) having 9 Si-C pairs (N = 9) in the unit cell and 
isosurfaces of charge densities of selected states at the T-point 
of BZ. Zero of energy is set at the Fermi level. Energy band 
gap is shaded. 



the edges. This results in a widening of the band gap. As 
for the other edge state band, it is located in the valance 
band. This band is also removed upon H-saturation, but 
the valence band edge is not affected. The band gap of 
H-saturated A-SiCNR(9) is direct. Energy bands and 
band decomposed charge densities are shown in Fig. [4] 

The variation of the energy band gap, Eq, of bare and 
H-passivated SiC armchair nanoribbons with the width 
of the ribbon N are calculated for 5<N<21 as presented 
in Fig. [5] For bare A-SiCNR's, the band gap is relatively 
smaller due to edge states as shown in Fig. [4] namely 
E G ~ 1.29 eV for N=5, but increases to ~ 2.4 eV for 
N=2\. Upon H-passivation of the dangling bonds at the 
edges, the edge states disappear and the band gap in- 
creases and gets direct. For N=5 the direct band gap is 
around 2.5 eV, but increases with N and eventually be- 
comes 2.38 eV for N=21. Hypothetically, the band gap is 
expected to reach the value of 2D planar SiC (Eg=2.53 
eV calculated within GGA). The main difference with 
graphene is that the band gap of bare armchair graphene 
nanoribbon decreases with increasing N and eventually 
vanish as N — > oo. Interestingly, H-saturated armchair 
SiCNRs with, N — 1, N and N + 1, exhibit a family 
behavior similar to one revealed in armchair graphene 
nanoribbonsPS For N > 15 the variation of Eq of H- 
saturated A-SiCNR is not significant. However, it should 
be noted that the band gap variation of both bare and H- 
saturated armchair SiC nanoribbons cannot be reconciled 
with the quantum confinement effects, since Eq increases 
with increasing N . This is due to other effects which over- 




Figure 5: (Color online) Optimized atomic structures for 
(a) Bare (b) H-passivated SiC nanoribbons. Large/blue and 
small/brown balls represent Si and C atoms respectively. The 
variations of band gaps of (c) Bare and (d) H-passivated SiC 
nanoribbons with the ribbon width N for 5 < N < 21. 



come the quantum confinement effect. The variation of 
band gaps of bare and hydrogen passivated armchair SiC 
ribbons are presented in Fig. [5] (c) and (d). In Ref. I3T1 
the variation of the band gap with N and resulting family 
behaviour is in agreement with present results. Whereas 
the same variation is not presented in Ref. [311 The zigzag 
SiCNRs exhibit interesting magnetic pro perties a nd were 
investigated extensively by earlier works BH 36 I 37 I For this 
reason, the zigzag SiCNRs are examined in the present 
study. 



VI. VACANCY DEFECTS AND ANTISITE 

It has been shown that the vacancy defects have re- 
markable effects on 2 D g raphene honeycomb structure 
and its nan oribbonspSHS Non-magnetic graphene sheets 
or nanoribbons can attain spin polarized states due to 
vacancy defects. We expect that similar effects of va- 
cancy defects can occur on the electronic and magnetic 
properties of SiC honeycomb structure. 



A. 2D Honeycomb SiC 

The effects of Si and C vacancies, Si+C-divacancy and 
C-Si-antisite are treated in periodically repeating super- 
cells. The size of supercell is optimized to allow negligible 
defect-defect interaction between adjacent cells. Here the 
width of the flat bands derived from the states of peri- 
odically repeating vacancies is taken as the measure of 
the strength of vacancy- vacancy coupling. A (7x7) su- 
percell is found to be suitable, since it leads to rather 
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Vacancies and Antisite in 7x7 2D SiC 



(a) C Vacancy (b) Si Vacancy 
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(c) Si+C Vacancy (d) Antisite 




H = ji = 

Figure 6: (Color online) Optimized atomic structure and mag- 
netic moment of vacancy defects calculated in a (7x7) super- 
cell of 2D SiC structure, (a) C vacancy; (b) Si vacancy; (c) 
C+Si divacancy; (d) C-Si antisite. In (b) the difference of 
spin up and spin down charges are shown. Large/blue and 
small/brown balls indicate Si and C atoms respectively. 



flat defect bands. At the same time it is not large and 
allows us to carry out numerical calculations within fea- 
sible computational time. The flat bands associated with 
vacancies can be considered as localized defect state (if it 
is in the band gap) or resonance state (if it is in the band 
continua). Here we consider localized defect states. Our 
results are presented in Fig. [6] for single C, Si vacancies, 
Si+C divacancy, C+Si antisite defects . 

A vacancy is generated first by removing a single atom, 
C or Si atom from each supercell of the monolayer of SiC 
as shown in Fig. [6] (a) and (b). Subsequently, the atomic 
structure is optimized. Single C vacancy in 2D monolayer 
of SiC is non-magnetic; Si atoms around vacancy with 
coordination number 2 are displaced in the transversal 
direction and do not induce any magnetic moment. 

As for Si vacancy, three C atoms around vacancy re- 
mained planar. Similar to the vacancies in graphene and 
BN, Si-vacancy induces a local magnetization in the sys- 
tem. Isovalue surfaces of the difference between spin up 
and spin down charge densities i.e. Ap^ clearly shows 
a spin polarization around the vacancy and a net mag- 
netic moment constructed therefrom. The calculated to- 
tal magnetic moment is 4 fig per supercell. The Si- 
vacancy in a repeating (7x7) give rise to defect states 
in the band gap. As for Si+C divacancy in Fig. [6] (c), it 
is again non-magnetic since the spins are paired; two C 
atoms around the vacancy choose to make a bond with 
each other. The band gap is also modified. Finally, we 



Table IV: Magnetic moments and positions of defect-induced 
state energies relative to the top of the valance band. Spin- 
up^) and spin-down(^) states are indicated. E (F) indicate 
whether the defect state is empty (full-occupied). 

Si+C Vacancy C- Vacancy Si- Vacancy Antisite 



M 








4 Ms 





Ej 


0.49 (F) 


0.15 (F) 


0.09 (F t) 


0.59 (F) 


E 2 


1.40 (E) 


1.82 (E) 


0.11 (F t) 


2.12 (E) 


Ez 


2.45 (E) 


1.95 (E) 


0.31 (F |) 




E± 






0.33 (F \) 




E 5 






0.38 (E I) 




E 6 






1.47 (E 4.) 




E-7 






1.49 (E I) 





consider the antisite defect. The resulting relaxed struc- 
ture is given in Fig. [6] (d). Lattice is distorted as C-C 
bond is shorter than Si-Si bond in the antisite case. It is 
noted that the calculated magnetic moments for single Si- 
and C- vacancy do not agree with Lieb's theorem™ which 
normally predicts lps net magnetic moment both for Si 
and C vacancies in Fig. [6] We attribute the discrepancy 
between the results of first principles calculations and 
Lieb's theorem to the structural relaxation occurred after 
the generation of vacancy and significant charge transfer 
from Si to C. The localized electronic states associated 
with the vacancy defects and antisite are deduced from 
the band structure calculations as presented in Table |TV| 
where the energies are given from the top of the valance 
band. 



B. Vacancy defects in SiC Nanoribbons 

The effects of the vacancy defects on the electronic and 
magnetic properties are treated for of H-passivated A- 
SiCNR(9) using a (4x1) repeating supercell. Our main 
motivation was to investigate what differences would oc- 
cur in a ribbon. Optimized structures, calculated total 
magnetic moments are presented in Fig. [7J Overall effects 
of vacancy, divacancy and antisite effects are similar to 
those in 2D SiC monolayer structure except in the anti- 
site case the exchanged Si atom moves out of the plane 
about 0.9 A. One localized state below the conduction 
band edge and one other state above the valence band 
edge occur due to antisite as donor and acceptor states, 
respectively. 



VII. FUNCTIONALIZATION OF SIC 
HONEYCOMB STRUCTURE BY ADATOMS 

Specific adatoms can bind to 2D SiC monolayer with 
significant binding energies. Adatom adsorption or dec- 
oration, as well as substitution for Si or C atoms in 
the honeycomb structure by foreign atoms can modify 
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Vacancy and Antisitc in A-SiC NR (9) 
(a) C Vacancy (b) Si Vacancy 




^ = 
(c) Si+C Divacancy 



M = 4 M B 
(d) Antisitc 



Table V: Eb, total magnetic moments fi, optimized heights 
of adatoms from the 2D monolayer of SiC, h and localized 
states of occurring in the band gap are given for each type of 
adatom. Empty (E), full (F), spin-up and spin-down states 
are indicated. Energies of adatom induced localized states are 
given relative to the top of the valence band. 






Co 


Fe 


P 


Ti 


Position 


HS 


HS 


BS 


TS 


E B (eV) 


2.099 


1.928 


1.758 


2.605 


M (Mb) 


1 


2 


1 


2 


h (A) 


1.392 


1.465 


1.743 


1.355 


Ei 


0.13 (F t) 


0.30 (F t) 


0.17 (F t) 


0.43 (F t) 


E 2 


0.23 (F t) 


0.31 (F t) 


0.16 (F t) 


0.60 (F |) 


Ez 


0.33 (F |) 


0.82 (F |) 


0.20 (F |) 


1.62 (F t) 


Ei 


0.41 (F |) 


0.82 (F |) 


0.20 (F t) 


1.62 (F t) 


E 5 


0.67 (F t) 


1.17 (F |) 


0.31 (F |) 


1.99 (E t) 


Ee 


0.80 (F |) 


1.93 (E t) 


0.33 (E |) 


2.28 (E t) 


E 7 


1.26 (E |) 


2.20 (E |) 




2.40 (E |) 



^ = ^ = 

Figure 7: (Color online) Magnetic moment of vacancy de- 
fects calculated in a (4 x 1) supercell of quasi ID SiC armchair 
nanoribbon with N=9, i.e. A-SiCNR(9). Dangling bonds at 
both edges are saturated by hydrogen atoms, (a) C vacancy; 
(b) Si vacancy; (c) Si+C divacancy; (d) C-Si antisite. In (b) 
the difference of spin up and spin down charges are shown 
while others are non-magnetic. Large/blue and small/brown 
balls represent Si and C atoms respectively. 



E B = E T [SiC + adatom]— E T [SiC (bare)]— E T [adatom]. 
We found that all of these adatoms are bound with a 
significant energy, which is larger than 1 eV. Flat bands 
indicate that states induced by the adatoms are rather lo- 
calized and hence adatom-adatom interactions are negli- 
gible. Therefore, adatoms treated here in supercell geom- 
etry can represent single (isolated) adatom. Optimized 
adsorption sites, total magnetic moments, heights from 
the SiC plane are given in Table |V| 



the properties of 2D SiC monolayer and its nanoribbons. 
This way, SiC honeycomb structures can be functional- 
ized. Adatom adsorption and substitution are considered 
within the periodically repeating (7x7) supercell geome- 
try to minimize the interaction between them. The cal- 
culations performed in a larger supercell (10x10) have 
also showed similar results which is an indication that 
the coupling between adjacent defects is negligible. 

A. Adatom adsorption 

As for adatom adsorption, we considered Al, Co, Fe, N, 
P, Ti by placing each of them on four different positions 
in the (7x7) monolayer of SiC and then by fully relaxing 
the whole system. The initial positions of the adsorption 
are on top of silicon atom (TS) , on top of carbon atom 
(TC), at the center of hexagon (HS), above the middle 
of the Si-C bond (BS). The distance between adatoms is 
12 A. Spin polarized calculations are carried out to deter- 
mine the binding structure and binding energy. Whether 
the adatoms are bound to the surface are examined by 
calculating the binding energies of these six different in- 
dividual atoms in terms of the calculated total energies as 



B. Substitution of Si and C by foreign atoms 

Here we examined the substitution of single Si or C 
atoms in 2D SiC honeycomb structure by various foreign 
atoms. Namely B and N substituting C atom; Al, As, 
Ga, P substituting Si atom. Similar to adatom calcu- 
lations, the substitution process is treated within peri- 
odically repeating (7x7) supercell. Because of periodic 
boundary condition the localized states appear as flat 
bands. We found that B, N, As and P atoms have a net 
magnetic moment fi= 1 fig . Ga and Al do not create any 
spin polarization. In Table |VI| the substitutional impu- 
rity states in the band gap of 2D SiC and resulting net 
magnetic moments are given. The substitutional foreign 
atoms, namely B, N, As, and P have either 3 or 5 va- 
lence electrons and hence destroy the spin pairing in the 
perfect honeycomb structure. The substitution of C and 
Si by these atoms gives rise to a net magnetic moment 
of 1 fig. As a result, 2D SiC can be magnetized without 
doping by transition metal elements. 
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Table VI: Energies of defects states occurring in the band gap 
are given for each type of substituted atom. Empty (E), full 
(F), spin-up and spin-down states are indicated. Energies are 
measured from the top of valence band. 



B N As P 




1111 


E 1 


0.23 (F, t) 1-47 (F, t) 0.75 (F, f) 1.18 (F, t) 


E 2 


0.69 (E, |) 2.37 (E, |) 1.18 (E, |) 1.77 (E, 4.) 


VIII. 


DISCUSSION AND CONCLUSIONS 



We present a study on ID SiC chains, 2D monolayer 
of SiC in honeycomb structure and its armchair nanorib- 
bons, 3D bulk SiC. We carried out stability analysis 
of those materials. Two dimensional monolayer of SiC 
is an ionic compound with charge transfer from silicon 
atoms to carbon atoms and has a hexagonal lattice form- 
ing a planar honeycomb structure. The calculation of 
phonon modes results in all positive frequencies in BZ 
and indicates the stability of the planar honeycomb struc- 
ture. This situation is in compliance with the previous 
works JESEH where the honeycomb structures of binary 
compounds of group IV elements or III-V compounds 
are stable in planar geometry if they have an element 
from the first row, such as B, C, N. The 2D SiC is a 
non-magnetic wide band gap semiconductor. However, it 
acquires net magnetic moment when a single Si- vacancy 
defect is created or Si and C atoms of the honeycomb 
structure are substituted by As, P, B and N. Single C- 



vacancy, Si+C divacancy and Si-C antisite defects does 
not give rise to any magnetic moment in the system. It 
is shown that 2D SiC can be functionalized through Si- 
vacancy and adatom adsorption or substitution of Si C 
by foreign atoms. 

Armchair SiC nanoribbons are found to be nonmag- 
netic semiconductors. We revealed the variation of 
band gap with the width of both bare and H-passivated 
nanoribbons. The variation of the band gap exhibits 
also family behavior. However, the band gap of arm- 
chair nanoribbons are smaller than the band gap of 2D 
SiC for small N. Therefore the confinement effect seen in 
narrow graphene armchair nanoribbons does not occur 
here. 

In conclusion, our state-of-the-art first-principles cal- 
culations demonstrate that 2D SiC monolayer is stable in 
honeycomb structure, and hence it has a strong chance 
to be synthesized in future. Above results indicate that 
bare and H-passivated SiC sheets and armchair nanorib- 
bons can present interesting properties which can be uti- 
lized in nanotechnology. Creating defects through va- 
cancies, doping through adatoms and substitutional im- 
purities can functionalize SiC honeycomb structure and 
hence add new magnetic and electronic properties. 
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